Clusters in a magnetic toy model for binary granular piles 
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Results on a generalized magnetically controlled ballistic deposition (MBD) model of granular 
piles are reported in order to search for the effect of "spin flip" probability q in building a granular 
pile. Two different regimes of spin cluster site distributions have been identified, a border line q c (f3J) 
where J is the interaction potential strength. 



INTRODUCTION 



The physics of granular matter has drawn a great deal 
of information from percolation theory ideas yj. Rela- 
tively simple models based on analogies have been imple- 
mented in order to describe realistic granular pile (static 
and dynamic) properties^, EH El IE El Bl ■ Segregation @, 
decompaction |9( and avalanches [10j point out to the 
existence of clusters. However, in describing such mate- 
rials it is crucial to consider that they are not made of 
symmetrical entities. It is necessary to introduce at least 
one degree of freedom for the grain with some coupling 
to an external field. One can imagine that the degree of 
freedom, is called a " spin" , coupled to a " magnetic field" , 
though the " spin" can represent any non-magnetic phys- 
ical feature of practical interest, like the grain roughness 
or shape feature. The spin role is to break the spatial 
isotropic symmetry. The direction of such a spin can rep- 
resent the position of a grain with respect to neighboring 
entities, as well as a rotation process. The grain-grain in- 
teraction can be imagined to be some (elastic-like) poten- 
tial containing information on the grain Young, rigidity, 
bulk modulus, and Poisson ratio, .. and geometric as- 
pects . Generalizations to more complex spin models 
are immediately imagined. 

Coniglio and Herrmann presented in |12| a related view 
of the granular packing problem and adapted the Ising 
model and Sherrington-Kirkpatrick spin glass model to 
granular phenomena obtaining two phase transitions in 
the system. The short range exchange energy J describ- 
ing a "spin-spin interaction" is analogous in granular 
matter to the contact energy between grains. An in- 
terpretation of J for flows can be found in Pandey et 
al.[13|]. A constrained Ising spin chain has also been re- 
cently considered and studied as a toy model for granular 
compaction (id). 

Something which is fully appreciated is the difference 
in constructing piles in presence or not of vertical walls. 
A pyramidal pile has not necessarily the same structure, 
arches, ... as a pile in a silo. Moreover it is quite unrealis- 
tic to build a pile from a single grain faucet. Finally, due 



to its anisotropic shape a grain can rotate during its fall 
e.g. in responding to the local wind, or difference in pres- 
sure between the upper and lower grain surfaces. Such 
remarks have motivated us into re-examining a magnetic 
Tetris-like model (or a magnetic rain) for which the grain 
is characterized by a " spin" which can flip during its fall 
under some energetic condition. 

In this paper the role of changing the depositing spin 
flipping probability during its fall in such a MBD (mag- 
netic ballistic depositon) model is shown to influence the 
pile density, the pile " magnetization" , and the cluster size 
distribution. In Sect EI we establish the algorithm rules 
and briefly comment upon them. In Scct lllll we present 
numerical results for the density and the magnetization 
of the pile. fSect HVlVjl . A critical percolation line is 
found separating two regimes for the size (mass) distri- 
bution of clusters. Finally, in Sect I VII a brief conclusion 
can be found. 



II. 



EXPERIMENTAL PROCEDURE 



The algorithm for the so called g-MBD model, in con- 
trast to the i-MBD model y~3, goes as follows: 

1. first, we choose a horizontal substrate of spins with 
a predetermined (for example antiferromagnetic- 
like or random) configuration; periodic boundary 
conditions are used; 

2. a falling (up or down) spin is dropped along one of 
the lattice columns from a height r max + 5a, where 
r max is the largest distance between an occupied 
cluster site and the substrate, 

3. at each step down the spin can flip i.e., change its 
"sign"; the "up" direction has a probability g; 

4. the spin goes down flipping until it reaches a site 
perimeter of the cluster at which time the local gain 
in the Ising energy 
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is calculated. The fall velocity is irrelevant and 
there is no backscattering. If the gain is negative 
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FIG. 1: The dependence of the density on q and f3J 



the spin sticks to the cluster immediately (sticking 
probability =1.0) and one goes back to step J5J). 
In the opposite case the spin sticks to the cluster 
with a rate exp(— A(3E) where Af3E is the local 
gain in the Ising energy. If the spin does not stick 
to the cluster it continues going down. Of course 
if the site just below the spin is occupied the spin 
immediately stops and sticks to the cluster. 



5. After the spin stops one goes back to step 0). 



In the i-MBD model[l5( a finite field H and q 
were assumed; in the present report we take H = 
enlarge the permissible values of q to [0, 1]. 



= 0.5 
but 



III. 



NUMERICAL RESULTS 



All results reported below are for a triangular lattice 
of horizontal size L = 100, when the pile made of clusters 
has reached a 100 lattice unit height, and after averaging 
over 1000 simulations. The substrate consists of spins 
with random direction. 



A. Density 

We define the density of a cluster as p = 

number of spins in the cluster . , . , n . , in 

r r-^n n — T-rr- — m wnicn obviously tne 

number of sites on tne lattice ■' 

number of lattice sites in the denominator = 10 000. 

FigG] illustrates the behavior of the density with re- 
spect to the q and (3J parameters. This figure con- 
vinces us that the results are symmetrical with respect 
to q = 0.5. In the i-MBD model the density varies be- 
tween 0.38 and 0.47 but in the present g-MBD a spread 
in density occurs, - from almost a completely compact 
pile in the AF case and large q to a loosely packed pile in 
the F case. The lowest density (0.38) occurs for (3 J —0 
and for border values of q, i.e. and 1, in the F-case. For 
strong enough positive interactions (/3 J > 4) the density 



FIG. 2: Dependence of the magnetization on q for several /3J 
values: F case (dotted lines), AF case (dashed lines) 



saturates toward the value p ss 0.47, in the F-case. In 
the AF case, the density varies between 0.38 and 1.0. 



B. Magnetization 

The dependence of the magnetization defined as 



M 



(2) 



is shown in Fig|3as a function of (3 J and q, where n + 
and ri- are the number of up and down spins respectively, 
i.e. n + = 10000/O+. M can be considered as a measure of 
the difference in grain orientations in the final packing. 

The magnetization appears to be the same in the re- 
gions (f3J >0,q< 0.5) and {(3 J <0,q> 0.5) (see Fig0. 
Notice that there is no terrace observed here in contrast 
to the i-MBD model case 



IV. 



PERCOLATION 



Define an up-spin percolating cluster as a cluster of 
up spins which extents from the bottom to the top of 
the sample. At fixed q and (3J the fraction p of piles 
consisting in such a percolating cluster of up spins was 
computed. The behavior of p with respect to q for several 
values of the (3 J parameter is shown in Fig|3J There is 
no percolating cluster for q < 0.75 in the F-case and 
0.85 in the AF-case. The q (qo((3J) ■= q when p > 0) 
dependence is shown in Fig^l Above go when p becomes 
finite a very fast growth of p is observed, at a given q c . 

The differences between AF and F are understood if 
one recalls that the density in the AF case is much more 
sensitive to a change in q than in the F case, for which 
a change of q induces only a very small variation of the 
density - recall that the density in this case is [0.38; 0.47]. 
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FIG. 3: Dependence of the percolation probability p c on as 
a function of q for different /3J values. 



For high q values the N s dependence stops to be an 
exponential law and becomes a power law (see log-log 
plots on FigEJ). Let us postulate 

N s (s)ocs~ kp for qpal. (4) 

In order to find the values of q for which the crossover 
effect occurs, i.e. from an exponential law regime (ELR) 
to a power law regime (PLR), we have estimated the 
slopes of the N s dependencies in both regimes. It seems 
that one should distinguish pile growth conditions with 
respect to q{fiJ) < q c and q((3J) > q c - 

Analogously to N s we can define Nh as the number of 
hole clusters. Unlike the spin clusters there is no differ- 
ence in Nh dependencies between low and high values of 
q, - all Nh, dependencies seem to follow a power law. 
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FIG. 4: Dependence of qo on (3J . 

On the other hand a similar variation in q values gener- 
ates piles with a wide interval of densities in the AF case. 

V. THE SIZE (MASS) DISTRIBUTION OF 
CLUSTERS 

The number of clusters with s size is called N s . Its 
behavior is shown in FigJSJfor low q values. Notice that 
the N s dependence is approximately a straight line on a 
semi- log plot (all logs are neperian), i.e. 

N s (s) cx e- kES for g«0, (3) 

where &e is a constant for a fixed (f3J,q) pair and s 
is the cluster size (mass). Observe that the existence of 
a large cluster is more probable in the F case than in 
the AF case. This is expected when one recalls that the 
F case favours neighbors with a similar spin direction 
during the deposition. 



VI. CONCLUSIONS 

We have presented the extension of the MBD model 
[l5| in order to find the role of the probability of spin 
flip during deposition (or the degree of freedom modifi- 
cation) in granular downward rain- like flow in 2D. The 
generalized model, hereby called g-MBD, is a nonequilib- 
rium ballistic deposition model with one degree of free- 
dom. One can imagine that the spins (grains) have differ- 
ent shapes and q can be related to a wind strength such 
that the grains favour one or another position in order 
to minimize the pile energy. We have examined the clus- 
ter properties through two order parameters, since two 
characteristic fields (J and q) are intrinsic to the model. 

We have investigated the size, or "mass", of the spin 
clusters created through simulation of the nonequilibrium 
deposition. The "quenching" of the degree of freedom on 
the cluster leads to two different regions of spin clus- 
ter geometric properties. In the low "field" q region the 
spin cluster mass distribution follows an exponential law, 
while in the high q region the distribution is character- 
ized by a power law. The transition between these two 
regimes is not sharp. The exponential law regime for 
high strength of interactions (i.e. \(3J\ > 4) seems to be 
universal, i.e. independent on the magnitude of the in- 
trinsic parameters but depends only on the sign of /3J, 
i.e. the characteristic contact interaction potential, - if 
they are ferromagnetic-like or antiferromagnetic-like, me- 
chanically repulsive or attractive. 
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semi log plot, low q, pj = 5.0 



semi log plot, low q, pj = -5.0 
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FIG. 5: Examples of the N s (s) dependence for low values of q on semi- log plots: (a) (3 J =5; (b) j3J =-5 



log-log plot, high q, pj = 5.0 



log-log plot, high q, pj = -5.0 
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FIG. 6: Examples of the N s (s) dependences for high values of the q parameter presented on log-log plots; (a) states for 
ferromagnetic-like case, i.e. f3J =5; (b) states for antiferromagnetic-like case, i.e. /3J —-5 
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